Dynamical suppression of 1// noise processes in qubit systems 
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We investigate the capability of dynamical decoupling techniques to reduce decoherence from 
a realistic environment generating 1/f noise. The predominance of low frequency modes in the 
noise profile allows for decoherence scenarios where relatively slow control rates suffice for a drastic 
improvement. However, the actual figure of merit is very sensitive to the details of the dynamics, 
with decoupling performance which may deteriorate for non-Gaussian noise and/or high frequency 
working points. Our results are promising for robust solid-state qubits and beyond. 

PACS numbers: 03.65.Yz, 03.67.Pp, 05.40.-a 



Noise processes characterized by a 1// power spectrum 
are ubiquitously encountered in Nature. While a unified 
theory of the underlying mechanisms remains elusive, 1/f 
noise plays a prominent role in dynamical phenomena as 
diverse as transport in electronic devices [l|, light emis- 
sion from astrophysical sources , statistics of DNA se- 
quences 0, and stock market prices In recent years, 
the continuous advances witnessed by device nanotech- 
nologies, along with the challenge to implement quantum 
information processing (QIP) in solid-state systems [3, 
have sharpened the demand for a detailed understanding 
of 1 // noise effects and for viable compensation schemes 
at the quantum level. In particular, 1/f noise due to 
fluctuating background charges (BCs) severely hampers 
the performance of single-electron tunneling devices 
and Josephson qubits in the charge regime [7|. 

Prompted by the experimental demonstration of a co- 
herent charge echo in a Cooper-pair box efforts are 
underway to explore the possibility of 1// noise reduction 
via active control techniques. Recent theoretical analyses 

^3 largely rely on deriving 1 // noise from the influ- 
ence of an harmonic oscillator bath e.g., through spin- 
boson models with sub-Ohmic damping [ll|. While this 
accurately represents environments consisting of a con- 
tinuum of weakly coupled modes, the inherent Gaussian 
distribution of the fluctuations fails at reproducing the 
distinctive properties oil/ f noise due to realistic discrete 
environments |l2l |. Compensation of non-Gaussian ran- 
dom telegraph noise (RTN) from a single bistable source 
is considered in |l3l |. However, an appropriate distri- 
bution of characteristic time scales is required to obtain 
genuine 1/f noise effects 0| , thus qualitatively changing 
the nature of the control problem. 

In this Letter, we present a comprehensive study of 
the effectiveness of decoupling methods 0, 0] at sup- 
pressing 1/f noise in a single qubit. Our approach fully 
captures both Gaussian and non-Gaussian effects for real- 
istic noise spectra and arbitrary operating points of the 
qubit. We find that the control performance depends 
critically on the frequency location of the dominating 1 // 
noise sources, the quality of the attainable suppression 



becoming comparatively higher as the latter shifts toward 
lower frequencies. For purely Gaussian 1/f dephasing, 
this implies the potential of significant coherence recov- 
ery (up to 75%) by using control rates which can be or- 
ders of magnitude slower than expected from the fastest 
characteristic time scale present in the system. 

Noise model.— A simple way for generating a 1// spec- 
trum is via an ensemble of M independent, randomly ac- 
tivated bistable processes. Let £,k{t) be an asymmetric 
RTN signal switching between values ±ffe/2 with rates 

7fc^\ Ik = Tfc^'' + 7fe~''- If a distribution P(7) oc I/7 is 
assumed for the switching rates 7fc € [7mjn, Imax] QllJI, 
the total fluctuation S(t) — J2k^h{t) exhibits a 1/f 
power spectrum of the form ^da;!) = A > 0, in a 

frequency range defined by effective cutoffs 7f„j„ ^ 7mm , 
Imax Imax- Wc focus On distributions of strengths Vk 
sufficiently peaked around their mean value {v) > 0, in 
which case A is proportional to the number Ud of fluctu- 
ators per noise decade, weighted by (w)^. 

The above phenomenological model applies to a vari- 
ety of 1/f noise processes e.g., due to delocalized charge 
traps or hopping defects in sohd-state devices Q|. In 
a fully quantum description, the RTN ensemble is re- 
placed by an environment E consisting of M two-state 
BCs 0, each coupled with strength Vk to the system 
and with strength Tk to an electronic band inducing re- 
laxation with rate 7^. Thus, He = X]fc ^k, with 
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Hk = ekblbk + n l^WkM + li.c] + 2^ eucliCki , (1) 
I 
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where euiSki are energy parameters and bk{b\), Cfcj(cJ,;) 
canonical fcrmion operators, respectively. The semi- 
classical approximation is accurate in a regime where 
each EC is strongly coupled to the corresponding band, 
Tk 3> Vk, implying the possibility to neglect quantum 
back-action effects from the system. Let gk = Vk/^k 
for each noise source. The resulting decoherence con- 
tribution is qualitatively different depending on whether 
gk ^ 1 (Gaussian source) or not [iJl- While purely 
Gaussian noise may be formally reproduced by an ap- 
propriate bosonic bath |lll |. non-Gaussian 1/f behavior 
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may exhibit non-equiUbrium saturation features and pro- 
nounced non-Lorentzian hneshapes |l2l| . 

Decoupling setting.— Our target system 5 is a single 
qubit, described by a Hamiltonian ifg = —['[laz+Aax]/'2, 
tr^-eigenstates {|0), |1)} defining the computational basis, 
and AE = •\/rP~f"A^ giving the energy scale of the free 
dynamics {h = 1). For the charge qubit realized in 0,01 j 
il is the separation between the two lowest charge states 
differing by a single Cooper pair, while A is the Joseph- 
son energy of the junction. Coupling with the BCs is in- 
troduced via an interaction Hse =X]fe ''^kb^bk <8) Cz, cor- 
responding semiclassically to the effective Hamiltonian 
HRTN{t) = S(t)(Tz. Both relaxation and dephasing dy- 
namics occurs, depending on the longitudinal or trans- 
verse nature of the fluctuations relative to the physical 
basis defined by the energy eigenstates. When A = 0, the 
physical and computational bases are aligned, and deco- 
herence is purely adiabatic. In the opposite case where 
51 = 0, purely non-adiabatic dissipation and dephasing 
take place, both mechanisms influencing decoherence in 
the z basis as the latter is now 90° displaced. Borrowing 
from the Cooper-pair-box terminology, we shall refer to 
these limiting situations as pure dephasing and charge 
degeneracy regimes, respectively [l^ . 

Dynamical decoupling strategies coherently average 
out the effects of unwanted interactions over a sufScie ntly 
long time scale by means of a tailored control field 
In the simplest formulation, decoupling is achieved by 
subjecting the system to cyclic sequences of instanta- 
neous (bang-bang) control pulses. We consider two el- 
ementary decoupling protocols specified by the following 
control cycles: Va = {Ai, tTj,, At, tt-^,} (asymmetric) and 
Vs = {Ai/2, TTj;, At, 7r„^, At/2} (symmetric), where At 
and ■K±x denote a free evolution interval of duration At 
and a controlled rotation of 90° about the ±x axis, re- 
spectively. While the two implementations are equiva- 
lent in the ideal limit of arbitrarily fast control At ^ 0, 
superior averaging is expected from the time-symmetric 
protocol in the realistic case of finite At thanks to the 
cancellation of higher-order corrections 0, 0, 0| . 

Decoupled dynamics.— Complete information about 
the dissipative qubit dynamics is contained in the ex- 
pectation values of the qubit observables ui, i = x,y,z, 
after evolution in the presence of the 1// disturbance 
with and without the decoupling field. For a generic Hs, 
the latter have been calculated numerically as the ensem- 
ble avera ge o ver the stochastic qubit state evolved under 
HRTN{t)m, {<7t{t)) =£{{'il^{t)\ai\iP{t))}, starting from 
known pure-state initial conditions = co|0)-|-ci|l), 

|coP + |cip = 1. We consider here two performance 
indicators: (cr+(t)) = {ax{t) + iay{t)) /2, starting from 
Co = ci = 1/ v^, which describes coherence in the z basis; 
and ((7z(t)), starting from cq = 1, ci = 0, which is a signa- 
ture for coherent oscillations in the x basis. At A = 0, an 
analytic benchmark has been obtained by both perform- 
ing the semiclassical approximation in the Heiscnbcrg 



equations of motion derived from Eq. |^ and by exactly 
solving the RTN dynamics. If Z{t) = \{a+{t))/ {a+{0))\ 
defines the controlled decay function, the following ex- 
pression holds in the presence of a single BC after a time 
t = 2N At corresponding to N control cycles : 
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Here, 2_Fi denotes the regularized hypcrgcometric func- 
tion 01, a = ^1 - g2 _ 2igtanh(e/2fcBT) 
z = 2e~'^^^ /a, and a, f{a) are respectively given by: 

p-iAt |- 

[l + g"^ + |ancosh(7a'At) 
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\a\'^ + {\ + g^ - 2ig{5p^q - 6pa)) 



+ 2{a' + 5poga") 
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\a\^ - [1 + g^ ~ 2tg{Speg - Spo)) 



Speq = tanh(e/2fcBr) and 5pQ — ±1 denoting the equi- 
librium and initial value of the BC population difference. 

For M noise sources, the individual contributions fac- 
torize, Z{t] = Y\t.Zk{t). In the absence of control, the 
result of (Eq. (4)) is recovered. In the continuous 
limit where At — t/N 0|, one can prove using 
Eq. (2) that .Zfe(t) 1 V/c, confirming perfect decoher- 
ence suppression for ideal decoupling. Thus, complete 
1// compensation requires access to control time scales 
S ^llmax, so as to queuch the influence of the fastest 
fluctuator present in the environment. We find that an 
almost full coherence recovery (better than 90%) can al- 
ways be achieved if decoupling is fast relative to the effec- 
tive upper cutoff, At < ^/jmax^ implying controls up to 
an order of magnitude slower. Tighter estimates of the 
minimum control rates able to ensure a significant de- 
coherence suppression are possible upon identifying the 
dominant noise sources in the relevant dynamical regime. 
While a detailed analysis is deferred to the salient 
features may be summarized as follows. 

Decoupling in the pure dephasing regime.— Adiabatic 
decoherence is insensitive to the energy scale AE — fl, 
but critically dependent on the coupling strength distri- 
bution, the overall time scale of the process being largely 
determined by (v) for small dispersions as assumed. An 
important consequence of AE being irrelevant is that de- 
coherence acceleration [l^ is never observed for A = 
with arbitrary control parameters. Symmetric decou- 
pling performs systematically better over its asymmetric 
counterpart, differences being larger for a number of con- 
trol cycles N < 5. For charge-echo sequences with realis- 
tic non-Gaussian noise and At values 0, the symmetric 
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readout signal can be up to 70% higher at times where 
the free coherence has entirely decayed (Fig. 1). Effects 
related to temperature as well as different initializations 
of the BCs (including sample to sample variations and 
correlations when repetitions of the readout process are 
needed as in |^) can be fully taken into account fioj . 
While this may be critical for a quantitative comparison 
with experiments, we anticipate no substantial changes 
as far as the decoupling efficacy is concerned. 
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FIG. 1: Symmetric {Vs, dotted) vs asymmetric {Va, dashed) 
decoupling for pure dephasing from a realistic 1// spectrum 
(Inset: jmin=l Hz (extrapolated), jmax = 10^^ Hz, Ud = 
1000, {v) = 9.2 ■ lO'' Hz, Av/{v) ~ 0.2). Upper panel: free 
evolution (solid line) vs echo signals. At = 1 ns. Lower panel: 
N — 5, At = 0.2 ns. Stroboscopic data points from Eq. 
@ are shown at t = 2NAt, whereas continuous curves result 
from averaging over 10"" RTN realizations. Each BC is initially 



assumed in a thermal mixture with e/2fc_BT ~ Spe 



0.08. 



A key factor determining the control effectiveness for 
dephasing processes is where (v) is positioned within 
[7mm,7max]- If noisc is purely Gaussian, Vk/jk < 1 Vfc, 
then (v) < "fmin for realistic spectra. Gaussian dephasing 
is dominated by the lowest spectral components, as wit- 
nessed by the diverging rate predicted by second-order 
perturbation theory for w — > 0. Hence, averaging of the 
slowest decades near 7mm suffices for a dramatic coher- 
ence improvement (Fig. 2(a)). Non-Gaussian behavior 
may arise either for a lower 7mm at fixed (v) (Fig. 2(b)), 
or for a higher (v) in a fixed range [jmimlmax] (Fig. 2(c)). 
Irrespective of whether saturation effects occur in the dy- 
namics of non-Gaussian charges p^ll9l |. it turns out that 
a substantial decoherence contribution originates from 
noise decades near (v), thus making their compensation 
essential. In practice, a necessary and sufficient crite- 
rion to obtain between 60-75% coherence recovery is to 
choose At < min(l/10«7™m , l/10'?(w)), with q = 1 or 2 
at most. For fixed ^mim this allows for pulse repetition 
times which can be longer and longer the lower the value 
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FIG. 2: 1// suppression under pure dephasing conditions, 
n = 1,A = 0. Top to bottom: (a) Purely Gaussian spec- 
trum, 7min = 10~*, -ymax = 100, Ud = 100, Iv) = 10~*; 
(b) and (c) Non-Gaussian spectra with parameters as in (a) 
except that 7min ~ 10^® in (b), and {v) — 0.01 in (c), respec- 
tively. 7maa; = 10 in all cascs. Solid lines depict free evolu- 
tions. Control parameters are: (a) and (b) At = 1000 (thin 
solid). At = 100 (dashed). At = 10 (dotted); (c) At = 10 
(thin solid). At = 1 (dashed). At = 0.1 ^ l/Tmax (dotted). 
Results are averages over 2 • 10* RTN realizations under Vs 
protocols. Each BC is initially in a pure state randomly sam- 
pled according to {5po) = 5peq — 0.08. Asymmetric RTN and 
analytic results from Eq. Q are also shown for N — 5 cycles. 



of (v) . In a strongly non-Gaussian regime where the lat- 
ter moves to higher frequencies, the advantages of slow 
decoupling are lost, and control rates ruled by Jmax^lmax 
may become necessary. 

Decoupling in the charge degeneracy regime.— A dis- 
tinctive feature of the dissipative dynamics in this limit is 
the sensitivity to the energy scale AE = A. For Gaussian 
noise, this may be understood from the fact that both re- 
laxation and dephasing arc governed (up to a factor 2) by 
a single rate, which depends only on the power spectrum 
S{u! = AE) to second order in the couplings. We find 
that the dynamical role of the energy scale AE extends 
beyond the Gaussian limit. In particular, irrespective of 
the Gaussian or non-Gaussian nature of the spectrum, 
decoherence rates at charge degeneracy do not substan- 
tially differ from adiabatic ones if the qubit operates at 
sufficiently low frequencies, AE — > "fmin- Shifting the 
working point to higher frequency has the effect of filter- 
ing out the majority of noise contributions, leading to co- 
herence times which are orders of magnitude longer than 
in the corresponding dephasing configuration (Fig. 3a). 
However, a hallmark of control in this regime is the pos- 
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FIG. 3: 1/f suppression under charge degeneracy conditions, 
n = 0, A = 1. Noise spectrum as in Fig. 2(c). Top panel, (a); 
coherence in the z basis for free (dashed) and controlled dy- 
namics. At = 10 (solid), At = 1 (dash dotted). The dephas- 
ing curve at AE = 1 is reproduced from Fig. 2(c) (dotted). 
Note the deteriorated performance for At = 10. Lower panels; 
damping (b, no control), acceleration (c, At = 10), and par- 
tial recovery (d, At — 1) of coherent charge oscillations. Com- 
plete recovery (not shown) is found for At = 0.1 ~ l/7maa;- 
RTN averages {Vs protocols) are taken over 2 ■ lO'' realiza- 
tions, with environment initialized as in Fig. 2. 

sibility to accelerate decoherence if cycle times long with 
respect to the free period 27r/Ai? are used. This is evi- 
dent in the oscillatory dynamics of {<Jz{t)) (Figs. 3b-d). 
If At ^ tt/AE, the improvement in the controlled ampli- 
tude does not compensate for the effective decay due to 
the frequency mismatch relative to the free oscillations, 
and overall acceleration results approaches 
tt/AE from above, the relative importance of these two 
mechanisms reverses, and a crossover to noise suppres- 
sion occurs at times which are shorter the shorter At. 
Full coherence recovery is guaranteed for At < 1/lOAi? 
which, however, may be close to l/'y^ax if large. 

Interestingly, a trade-off emerges in the charge degen- 
eracy regime between noise effects which are stronger but 
easier to decouple (as AE shifts toward low frequencies); 
and noise effects which are substantially reduced from 
the beginning but harder to suppress (as AE increases). 
A similar conclusion is likely to hold for generic qubit 
parameters. Thus, while it may seem counterproductive 
to consider low operation frequencies in the presence of 
1/f noise, in practice both the relevant noise level and 
the available control resources should guide the choice of 
a working point able to maximize the control efficiency. 



Conclusion.— We have characterized the performance 
of decoupling techniques at reducing 1/f noise based on a 
realistic model. Beside identifying noise scenarios where 
decoupling is highly effective with affordable rates, our 
results suggest the possibility of using control as a di- 
agnostic tool to infer spectral properties (such as (v) or 
lmax)^ which are not directly measurable. While actual 
details will be device-dependent and require a dedicated 
analysis incorporating the underlying physics, we expect 
that our main conclusions will have a wide range of appli- 
cability, including noise spectra with power-law behavior 
of the form l/|aj|" with a > 1 0|, or alternative qubit de- 
sign as demonstrated in . Ultimately, our results hold 
the promise that decoupling methods may substantially 
improve the prospects for reliable solid-state QIP. 
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